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Commuting ordinary differential operators with polynomial 
coefficients and automorphisms of the first Weyl algebra 
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Abstract 

In this paper we study rank two commuting ordinary differential operators with polyno¬ 
mial coefficients and the orbit space of the automorphisms group of the first Weyl algebra 
on such operators. We prove that for arbitrary fixed spectral curve of genus one the space of 
orbits is infinite. Moreover, we prove in this case that for for any n > 1 there is a pair of self- 
adjoint commuting ordinary differential operators of rank two L4 = (df + V(x)) 2 + W(x) , 
Lq , where W(x), V(x) are polynomials of degree n and n + 2 . We also prove that there 
are hyperelliptic spectral curves with the infinite spaces of orbits. 


1 Introduction 

The group of automorphisms of the first Weyl algebra A\ = {^” =0 Uj{x)d 3 x , Uj G C[x]} acts on 
the set of solutions of the equation 

n 

f(X, Y)=J 2 = 0 , A,, a l3 G C, ( 1 ) 

j,i=0 

i.e. if X,Y G A\ satisfy (pQ) and cp G Aut(A \), then <p{X),ip{Y) also satisfy flU). The group 
Aut(A\) is generated by the following automorphisms 

<Pi{x) = ax + f 3 d x , <pi(d x ) = 72; + Sd x , a,7, <5 G C, ad - ^7 = 1 , 

<p 2 (x) = x + Pi(d x ), P2{d x ) = d x , 
tp 3 (x) = X, ip 2 (d x ) = d x + p 2 (x), 

where Pi, P2 are arbitrary polynomials (see m )• So, Aut(Ai) consists of tame automorphisms. 
A natural and important problem is to describe the orbit space of the group action of Aut(Ai) in 
the set of solutions of (|TJ) . If one describes the orbit space it gives a chance to compare End(A\ ) 
and Aut(Ai) (End(Ai) consists of endomorphisms ip : A\ — > A\ , i.e. [<p(d x ), <p{x)\ = 1 ). 
Let us recall the Dixmier conjecture: End{A\) = Aut(Ai) , or in other words, if differential 
operators L n ,L m with polynomial coefficients satisfy the string equation 

[Lru L m ] — 1, 

then L m ,L n can be obtained from x,d x with the help of compositions p>j above (the general 
Dixmier conjecture for A n is stably equivalent to the Jacobian conjecture due to ED- Berest 
has proposed the following interesting conjecture: 

If the Riemann surface corresponding to the equation f = 0 with generic otij G C has genus 
g = 1 then the orbit space is infinite, and if g > 1 then there are only finite number of orbits. 

*The first author was supported by RSF (grant 14-11-00441). 
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One can prove that if there are finite number of orbits for some equation JT]) then End(Ai) = 
Aut(Ai) . 

In this paper we consider the equation 

Y 2 = X 2a+1 + c 2g X 29 + ■ ■ • + ciX + c 0 , X, Y € Ai , Cj e C. (2) 

Using Schur’s arguments [ 3 ] one can prove that if X. Y £ A\ satisfy (J 2 J) then XY = YX . 
Our approach to the above problem is based on the Krichever-Novikov theory of commuting 
higher rank ordinary differential operators. Let us recall some basic notions and facts related to 
commuting differential operators. If L n = )UJ =0 v 3 i x )^i > = J2T=o Uk ( x )^x commute then 

there is a Burchnall-Chaundy’s polynomial F(z, w) which vanishes the operators, F(L n , L m ) = 
0. 

The spectral curve T defined by the equation F = 0 is irreducible and is completed at 
infinity with a unique point q . The spectral curve parametrizes common eigenvalues of L n and 
L m , i.e. if L n ip = zip, L m ip = imp, then (z, w) £ T. The dimension of the space of common 
eigenfunctions for generic P = ( z,w ) € T is called the rank. Commutative rings of ordinary 
differential operators were classified by Krichever a, a- In the case of rank one eigenfunctions 
are Baker-Akhiezer functions, found by Krichever. The case of rank l > 1 is very complicated. 
In this case the eigenfunctions can not be found explicitly. Operators of rank two corresponding 
to elliptic spectral curves were found by Krichever and Novikov [6], operators of fourth order 
have the form 

Lkn = (d 2 x + u) 2 + 2c x (p(7 2 ) - p(m))d x + (c a ,(p( 7 2 ) - p(7i)))* - £>( 72 ) - p(7i). 


where 71 (x) = 70 + c(x), 7 2 (z) = 70 - c(x), 
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<x) = -jp I + 1 ^f + 2d>(7i,72)c, - ^ 70 + c, 7 o - c) - ^(71,72)), 


$(71,72) = C(72 -71) + C(7i) - c(72), 

C{z),p(z) are the Weierstrass functions, c(x) is an arbitrary smooth function, 70 is a constant. 
The operator Lkn commutes with a six order differential operator Lkn ■ 

Let us formulate our main results. 


Theorem 1 . 1 . For arbitrary integer m > 0 and arbitrary spectral curve T given by the equation 
w 2 = z 3 + c 2 z 2 + c\z + Co there are polynomials 

= o m _|_ 2 x T • • • T olq, YL m = !3 m x + ... + (3q, ct m _|_ 2 ^ 0 , (3 m ^ 0 

such that the operator 

L4,m = (d 2 + V m (x)) 2 + W m (x) 

commutes with a six order operator L^ m . The spectral curve of L^ m ,L§^ m coincides with T . 
At m = 1 we have 

L41 = ( d 2 + a 3X 3 + a 2 x 2 + a\x + ao) 2 + 20307 03 7^ 0. 

At 03 = l,oi = a 2 = 0 the operators coincide with the Dixmier operators jT]. The 

example of Dixmier was the first example of commutative subalgebra in A\ . It is an interesting 
problem how to obtain L4 :Tn ,LQ trn from Lkn, Lkn 2 At m = 1 the answer is given in the 
Grinevich’s theorem [ 7 ]: 
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• Operator Lkn corresponding to the curve w 2 = 4 z 3 + g 2 z + g 3 has rational coefficients 
if and only if 

/ \ _ r dt 

C[X) U)y/*t* + 92 t + 9 s’ 

where q(x) is a rational function. If 70 = 0 and q(x) = x, then Lkn coincides with 
£4,1 • 

Theorem 1.1 allows to prove the following theorem. 

Theorem 1.2. The set of orbits of the group Aut(A \) in the space of solutions of arbitrary 
equation 

Y 2 = X 3 + c 2 X 2 + Cl X + c 0 , X,7e A uCj £ C 

is infinite. 

Commuting operators of rank two of order 4 and 4 g + 2 corresponding to hyperelliptic 
spectral curves of genus g were studied in [8]. With the help of methods of [8] one can construct 
rank 2 operators at g > 1 . For example 

l\ = (d 2 + a 3 x 3 + a 2 x 2 + a\x + a 0 ) 2 + 3(9 + 4 )a 3 x, a 3 / 0 

commutes with an operator l \ 9+2 [H] - Mokhov [ 9 ] proved that if one apply elements of Aut(Ai) 
to l\, L; g+ 2 then one can obtains operators of rank l = 2 k and l = 3 k , where A: is a positive 
integer. For example if we apply the automorphism <p(x ) = d x ,cp(d x ) = —x to l\,l \ 9+2 we 
obtain rank 3 operators. Herewith 

ip(L[) = (a 3 d 3 + a 2 d 2 + aid x + a 0 + x 2 ) 2 + g{g + 1 )a 3 d x . 

Another important example constructed in m is the following. The operator 

l\ = (d 2 + ai cosh x + a 0 ) 2 + at\g(g + 1) cosh x, a\ / 0 

commutes with L \ g+2 . Using l\, L \ g+2 Mokhov constructed examples of operators of arbitrary 

rank l > 1 m (we discuss this construction in section 2 ). Let be a spectral curve of l\, l \ 9+2 
given by the equation 

w 2 = z 2g+1 + c\ g z 29 H-1- c\z + Cq. ( 3 ) 

Coefficients Cj can be found with the help of a recurrent formula (see Lemma 1 bellow). Probably 
for all g the curve is not singular for general set of parameters ao,ai . For small g using 
Lemma 1 one can check this by direct calculation. 

Theorem 1.3. The set of orbits of the group Aut{A\) in the space of solutions of the equation 
Y 2 = X 2g+1 + c\ g X 2a +-1- c\x + cj, X,Y £ A\ 

is infinite. 

It would be interesting to check the Berest conjecture at g > 1 for generic equation (HI) 
having a nonconstant solution in A\ . 

Remark 1.1. The group Aut(A\ ) acts on the set of rings of commuting differential operators 
with affine spectral curves considered in Theorems 1.2 and 1 . 3 . One can prove that the space of 
orbits is also infinite. 
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2 Method of deformation of Tyurin parameters 


Every ring A of commuting ordinary differential operators is isomorphic to a ring of meromor- 
phic functions on spectral curve F with a pole in some point q € T (we consider in this section 
the case when T is nonsingular, i.e. T is a Riemann surface). For a meromorphic function 
f(P ), P G T with pole in q of order n we have Lfi/)(x,P) = f(P)i/j(x, P) where Lf e A is 
a differential operator of order In , l is the rank of commuting operators, if) = (ipi, ..., ipi) is a 
vector Baker-Akhiezer function. Function t/j can be reconstructed from the following spectral 
data (see @D 

{F,q,k x , 71, ..., 7 Zg , ai, ..., cui(x), ..., (x)}. 

Hire k~ l is a local parameter near q , g is the genus of T, Xj € T , atj = (a^i,..., 07,2-1) 
is a vector, uij(x) is a smooth function. The set (7, a) is called the Tyurin parameters. This 
parameters define a semi-stable holomorphic rank l vector bundle on T of degree Ig with 
holomorphic sections 771, ..., 7/ . The points 71, ..., 7 i g are points of their linear dependence of 
the sections 

l -1 

i= 1 

The vector-function ip is defined by the following properties. 

1 . In the neighbourhood of q it has the form 


ip(x,P) 


^ 2 ts{x)k 


vs=0 



<F(x,k), 


where £0 = (1,0,..., 0), £j(x) = ..., ^(x)), the matrix satisfies the equation 




( ° 

1 

0 

0 

0 

\ 



0 

0 

1 

0 

0 


A<F, 

A = 

0 

0 

0 

0 

1 




\ k + UJ 1 

0J2 

w 3 

• • • Ul- 1 

0 

/ 


2 . The components of -0 are meromorphic functions on T\{q'} with the simple poles Ji, ■ ■ ■ ,Ji g , 
and 

Res 7i 0 j = aijReSy^i, 1 < i < Ig, 1 < j < l — 1 . 

The main difficulty to construct operators of rank l > 1 is the fact that the Baker-Akhiezer 
function is not found explicitly. In the recent paper [ 12 ] were shown that the class of Baker- 
Akhiezer functions contains some known special functions. 

Let us recall the method of deformation of Tyurin parameters [6]. The main idea of this 
method is to study the linear differential operator which vanishes the common eigenfunctions. 
The common eigenfunctions of commuting differential operators of rank l satisfy the linear 
differential equation of order l 

^ l \x,P) = xo(x,P)ifi(x,P) -I-hXi-i(x,P)^ (Z_1) (x,P). 

The coefficients Xi are rational functions on F with the simple poles Pi(x), ..., Pi g (x) £ T , 
and with the following expansions in the neighbourhood of q 

Xo(x,P) = k + g 0 (x) + 0 (k~ 1 ), Xj{x, P) = gj{x) + 0(fe _1 ), 0 < j < l - 1, 

Xi-i(x,P) = 0(k~ 1 ). 
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Let k — 7 i(x) be a local parameter near Pi(x ). Then 


= + d i :j (x ) + 0 (fc - 7i(s)). 

Functions Cij(x),dij(x) satisfy the following equations (see [1]). 

<H,i-i(x) = -7,'(®), (4) 

d*,o0*0 = «i,o(a:)ai,!-2(3;) + «i,o 0*0*0,z-i 0*0 - a' ifi (x), (5) 

di,j{x) = aij(x)a ij i- 2 (x) - aij-i(x) + a i j(x)d ij i- 1 (x) - a' iJ (x),j > 1, (6) 

where aij(x) = c c (’ J , 0 < j < l — 1, 1 < i < Ig. To find Xi one should solve the equations 
([4])- (( 6 |). Using Xi one can find coefficients of the operators. At g = 1, 1 = 2 Krichever and 
Novikov [5] solved these equations and found the operators L^n ■ Operators of Krichever- 
Novikov and it applications were studied in m-m Operators of rank 3 corresponding to 
elliptic spectral curves were found by Mokhov [20] , In JT) [271 some examples of operators of 
rank 2,3 corresponding to spectral curves of genus 2-4 were constructed. 

In® commuting operators of rank two of order 4 and 4g + 2 corresponding to hyperelliptic 
spectral curves were studied 

Lx 0 = zip, L 4 g+2 ip = wip, w 2 = F g (z ) = z 29+1 + C 2 9 z 29 4-h c 0 . 

Common eigenfunctions of L 4 and L 4 g+2 satisfy the second order differential equation 
ip" - Xi(*. P W - Xo{x, P)ip = 0, P = (z, w ) € r, 
where Xo(x, P), Xi(x, P) are rational functions on T satisfying equations ([1]) (l 6 |) . 

Theorem 2 (®) The operator L 4 is formally self-adjoint if and only if 

Xi (x,P) = Xi(x,cr(P)), 
where a is the hyperelliptic involution on T . 


Theorem 3 (®) If I /4 is formally self-adjoint, i.e. L 4 = (d% + V(x)) 2 + W(x), then 

1 Qxx W jr Qx 

X °~~2 ~Q~ + Q~ ’ Xl ^-Q’ 

where Q = z 9 + a g -i(x)z 9 ~ 1 + ■ ■ ■ + ao(x), ao(x),..., a g -\{x) are some functions. The function 
Q satisfies the equation 


4 Fg{z) = 4 (z - W)Q 2 - 4 V(Q X ) 2 + ( Q xx ) 2 - 2 Q X Q XXX + 2Q(2V X Q X + 4 VQ XX + d*Q). (7) 


From Theorem 3 it follows 

Corollary 1 The function Q satisfies the linear equation 

d x Q + 4U Q xx x + 6V x Qxx + 2(22 — 2IU + V XX )Q X — 2 W X Q = 0. (8) 
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Corollary 2 If g = 1 then 


V = 


-16Fi(±(-c2 - IT)) + W 2 * - 2 W X W X 
4W? 


(9) 


where F\ defines the spectral curve w 2 = F\ (x) = z 3 + C 2 Z 2 + c\z + cq. 


With the help of Theorem 3 many examples of rank 2 operators were recently constructed 
(see [25j [27])- 

Let us consider commuting operators l\ g+2 [TO]. These operators do not commute with 
operators of odd orders [28j, hence these operators are operators of true rank 2. The polynomial 
Q for L\,l\ g+2 has the form (see [TO]) 


Q(x, z) = A g (z ) cosh s x + • • • + Ai(z) cosh x + Aq(z), 


where 
A* = 


1 


-&4 


8(2s + l)ai(g(g + 1) - s(s + 1)) 

( s + 2)! 


4A, 4 . s L±2! - 8A, i .,h±Sl(2 ao + + 4s + 5)- 


S! 


s+2" 


S! 


(2s + 3)ai + 4A s _|_i(s + l)((s + 1) (4ag + (s + 1) + 4 z) I , 0 < s < g, (10) 


we assume that A s = 0 at s < 0 and s > g , A g is a constant. 

Lemma 1 ([8]) The spectral curve T^ of L\,L\ g+2 is given by the equation 

w 2 = F g (z) = — (4 AqZ — 4^4o^4i«i — 16j4 2 (ao + 1) + 48^4) + 4ao^4i + 4A 2 — 2A±(6A^ — ^4i)) , 
where Aj(z) are defined in m>- 


Examples: 

!) 9 = 1 

F 1 (z) = z 3 + (- - 2 a 0 )z 2 + — (1 - 8a 0 + 16ag - 16 a\)z + 

z lo 4 

2 ) g = 2 , let for simplicity of formulas ao = 0 

F 2 (z) = z 5 + — z 4 + — (321 -336a?)z 3 + ^-(34-531a?)z 2 + (1 - 189a? + 108a? )z + 24a? + 513af. 
2 16 4 

The spectral curves defined by the above equations are not singular for the general parameters. 
Mokhov mi found a remarkable change of variable 

x = In (y + y/y 2 - l) r , r = ±l,±2 ,..., 

which reduces the operators L\,L\ g+2 to the operators with polynomial coefficients. In partic¬ 
ular, l\ in new variable y gets the form 

L\ = ((1 - y 2 )dy - 3 yd y + aT r (y ) + b) 2 - ar 2 g{g + l)T r (y), afi 0, 

b is arbitrary constant, T r (y ) is the Chebyshev polynomial of degree |r| . Recall that 

T 0 (y) = 1, Ti(y) = y, T r (y) = 2yT r _i(y) - T r _ 2 (y), T_ r (y) = T r (y). 
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Chebyshev polynomials are commuting polynomials, i.e. 

Tn{T m (y)) = T m {T n {y)) = T n _|_ m (y). 

If one applies the automorphism 

<p(y) = ~9 y , tp(d y ) = y, if <E Aut(Ai) 

to the operators l\, L \ g+ 2 written in y variable, then one gets operators of orders 2r, (2 g + 1 )r 
of rank r nn and 

<p{L\) = ( aT r (d y ) - y 2 dy - 3 yd y + y 2 + bf - arg(g + 1 )T r (d y ). 

3 Proof of Theorems 1.1—1.3 

3.1 Proof of theorem ll.il 

Let us rewrite ([9|) in the form 


m 2 V = -16Fi(^(c 2 - W)) + wl x - 2W X W XXX , (11) 

Note that from (fTTh it follows 

- 4 FI(±(-c 2 - W)) + 2 V X W X + 4VW XX + W xxxx = 0. (12) 

Further we assume that V, W are polynomials 

\ = OL n x n + ... + ao, W = /3 m x m + ... + /3o, 0, Pm 0- (13) 

Equation (flTT) is equivalent to the system of equations: equation (fT2l) and the equation on free 
terms of m which is 

a 0 pl = - fa)) + Pi - 3/?i/? 3 . (14) 

Let us prove the following important proposition. 

Proposition 3.1. For any m > 0 there exists a solution of the equation dill) of the form H13\) . 
where n = m + 2 . 

Proof. Equation d!2D is equivalent to a system of 2 m + 1 equations in 2m+ 4 variables a^Pj . 
Note that all equations have degree 2 and the set of their solutions consists of points in C 2m+4 
(with coordinates a*, Pj ) which lie in the intersection of 2m +1 quadrics defined by these equa¬ 
tions. By [29^ Ch.l,Th.7.2] the intersection X of these quadrics in P 2m+4 (with homogeneous 
coordinates cti,Pj,u) is non-empty and each its irreducible component has dimension greater 
or equal to 3. By the same reason the intersection of X with the hyperplane Z = {u = 0} at 
infinity is non-empty and each its irreducible component has dimension greater or equal to 2. 

To prove the proposition it is sufficient to prove that for any fixed m > 0 there is a two- 
dimensional irreducible component of In Z. From this fact we can conclude that affine part 
of the intersection of quadrics is non-empty. 

The homogeneous parts of our equations in P 2m+4 not depending on u can be easily written: 
these are exactly the coefficients at x l of the sum 

4VW XX + 2V X W X -31T 2 . (15) 
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Let us introduce the following notations: 

m+n—2 m+n—2 m+n—2 

i4w* = bixi > = E ^ ^ 2= E 

2=0 2=0 2=0 

Then the intersection TflZ is given by the equations 

4cj + 2 bi — 3 di = 0, i = 0,..., 2m. (16) 

Note that the coefficients 5j,Cj,dj can be written in the following form: 

2 2 2 

di ^ ( /3j—fc/3fc, hi ^ ' Bk,jC%i—k+lfik+1 ; Q ^ ' 6fcyOy_ kfik+2 ; (17) 

/c=0 fc=0 fc=0 

where = (A: + 1)(* — k + 1 ),Ck,i = (k + l)(/c + 2) are positive integers, and we set ay = 0 
if j > n , /3j = 0 if j > m. 

The next observation is: equations (|16lh (|17l) always have a solution of the form 

P = (a n + 0 : p m ± 0 : 0 : ... : 0) 

for any m > 0. Indeed, if ao = ... = a n -i = + = ... = (3 m -i = 0, then only 2m-th equa¬ 
tion from (1161) remains to be non-trivial, and this equation becomes a quadratic homogeneous 
equation linear in a n and quadratic in /3 m : 

(233m—1,2m T 4C m _2 J 2m,)o : n/3m 3/3^ — 0. 

Thus, we can set /3 m = 1 where from a n = 3/(2B m _i t 2m + 46 m -2,2m) ■ 

Let us prove that for any fixed m > 0 any irreducible component of X D Z containing P 
has dimension 2. 

If m = 1 then there are only 3 equations in (1161) 

46*0,00:0^2 + 2-Bo,oai/32 — 3/3 q = 0, 

4(Co,iai/3 2 + 6 * 1 , 20 : 0 / 33 ) + 2(5o,ia2^i + Si,iai/3 2 ) — 6/3o/3i = 0, 

4 ( 6 o, 202/?2 + 6 * 1 , 201^3 + 62 , 200 / 34 ) + 2(5o, 20 : 3^1 + -Bl, 202/32 + 732 , 201 / 33 ) ~ 3 ( 2 /?o ^2 + /3i ) = 0 , 
and their Jacobi matrix at P has the following form: 

/ * 2i3o,o/3i 0 0 * * \ 

I * * 22?o,i/3i 0 * * I ) 

\ * * * 273 0 ,2/3i * * J 

where the first columns denote derivations with respect to oo,..., 03 , and the last two columns 
denote derivations with respect to /3q, (3\ . The rank of the matrix is 3, so, these equations define 
a smooth variety in the neighbourhood of the point P of dimension two. 

For generic m the point P might not be regular. Nevertheless, any irreducible component 
containing P has a dense subset of smooth points. At any such point Q the Jacobi matrix J 
can be written in the following form. It can be divided in two blocks: one consists of m + 1 
columns (derivations of equations with respect to /3o,, /3 m ), and another one consists of n +1 
columns (derivations of equations with respect to a n ,... ,«o ). We shall describe only essential 
columns for us. 

The columns of the first block are (to save the space we shall write them as rows): 


2 -nd column: {jo t0 a 1 , (jo, 1^2 - 6/3 0 ), (jo, 20:3 - 6/3i),..., (jo, n _ia: n - 6/3 m ), 0 ,..., 0 ) 



3 -d column: (ji,oa 0 , Ji.iOi, (71,20:2 - 6 / 3 0 ),..., (ji, n otn ~ 6 / 3 m ), 0 ,... , 0 ) 

4 -th column: ( 0 ,72,100,72,201, (72,302 - 6 / 3 0 ),..., (72,71+1 o n - 6 / 3 m ), 0 ,... , 0 ) 

5 -th column: ( 0 ,0,73,20:0, .73,301, (73,402 - 6 / 3 0 ), • • •, (73,n+2a„ - 6 / 3 m ), 0 ,... , 0 ) 

(HI —|— 1)-1h column. (0, . . . , 0,7'm—l,m—200,7'm—l,m—10:1, (ira-l,m®2 6/?o), • • • , (im-l,2m®n 6 / 3 m)), 
the columns of the second block are: 

1-st column: (0, . . . , O, 7 o, m+ i/ 3 i, 7 i, m+2 / 02 , • • • ,jm-l, 2 mPm) 

2- nd column: (0, . . . , 0 ,j O ,mPl,jl,m+lP 2 , • • • ,jm-l,2m-lPm,0) 

3 - d column: ( 0 , . . . , 0 , jo,m-i/ 01 , 7 !,m/ 02 , ■ ■ ■ , 7 m-l, 2 m- 2 / 0 m, 0 , o) 


n-th column: (7 o,o/3i,7i,i/3 2 , • • • ,jm-l,m-lPm, 0,..., 0) 
(n+l)-th column: (ji, 0 / 3 2 , • • • ,7m-i,m-2/0m, 0,... ,0), 
where the numbers jkj are defined as 


3k,l — 2 + 4C'fc_ 1 ,i, 

where we assume B^i = 0 if k > l and C^i = 0 if k < 0 . 

Without loss of generality we can assume that the point Q belongs to a sufficiently small 
neighbourhood of the point P (in the complex topology), such that, for fixed numbers jk,i , 
the modules of all terms of the matrix J , except the terms containing /3 m = 1 and a n , are 
comparable with some 0 < e <C 1 (i.e. they are < e but > e 2 ). We call such terms comparable 
with e. 

We have the following possibilities now. If there is a smooth point Q such that its coordinate 
«o 7 ^ 0 or a\ 7 ^ 0 , then the rank of the matrix J is 2m+ 1 , i.e. the dimension of the component 
is two. Indeed, we can first apply the Gauss elimination algorithm to kill all terms of the right 
part of the matrix lying over terms containing /3 m . We can choose e small enough such that the 
terms of the left part of J will change, but the top non-zero elements of the first m — 2 rows 
will remain non-zero and comparable with e, and all elements over them will be comparable 
with e 2 . Applying again the Gauss elimination algorithm we can kill all elements in the columns 
except these top non-zero elements, thus obtaining 2m + 1 linearly independent rows in the 
matrix J. 

Note that the case ao = 0 , ai 7 ^ 0 (i.e. ao = 0 for all smooth points) is in fact impossible: 
in this case the whole component belongs to the hyperplane ao = 0 . But then the dimension of 
the component must be 1 , a contradiction. 

Now we claim that there exists a smooth point such that ao 7 ^ 0 or a\ 7 ^ 0 . Indeed, if there 
are no such smooth points, then the whole component belongs to the intersection of hyperplanes 
ao = ai = 0 (cf. 129 . Ch.l,ex.l. 6 ]). Note that in this case from 0 -th equation in (|16D it follows 
Po = 0 , and from the 1 -st equation it follows 02 /3i = 0 . 

Let’s show first that a 2 = Pi = 0. If there is a smooth point in the component with 
a 2 7 ^ 0 , then the Jacobi matrix of our system restricted to the 2m-dimensional intersection of 
hyperplanes ao = ai = Po = 0 reduces to the following matrix. 

The columns of the first block are (to save the space we will again write them as rows): 

1- st column: (jo,ia 2 , ( 70 , 20:3 - 6/3i),..., (jo,n-i&n ~ 6/3 m ), 0,..., 0) 

2- nd column: (0,7i, 2 a 2 , ( 71,303 - 6/3i),..., ( 7 i, n a n - 6/3 m ), 0,... , 0) 
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3-d column: (0, 0, .72,30:2, • • •, (j2,n+i«n - 6/3 m ), 0,..., 0) 


m-th column: (0 ,... , 0, j m -i, m a 2l ■ ■ •, ~ 6 P m ))\ 

the columns of the second block are: 

1 -St Column: ( 0 , . . . , 0 , jo,m+lPl, Jl,m+2/3 2 ) • • • , jm— 1 ,2m firm) 

2 - nd column: (0, . . . , 0, jo,mPl, jl,m+lP 2 , ■ ■ • , jm-l, 2 m-lPm, 0) 

3 - d column: (0, . . . , 0, jo,m-lPl, jl,mP 2 , ■ ■ ■ ,jm-l, 2 m- 2 Pm, 0,0) 


(n-l)-th column: (jo,iPi, ■ ■ ■ ,jm-i,mPm, 0 , • • • , 0 ), 

where the m columns of the first block denote derivations with respect to Pi,..., j3 m , and the 
n — 1 columns of the second block denote derivations with respect to a n ,... ,a 2 ■ 

Since f3\ must be equal to zero, we can apply the same arguments as above and obtain that 
the rank of this matrix is 2m . But this is impossible, because the dimension of the component 
is not less than two. 

If the whole component belongs to the intersection Y = {ao = at\ = a 2 = Po = 0} , but there 
are smooth points with /3i / 0 , then the 2 -th equation in (fl 6 l) reduces to 2 i?o, 2 « 3 / 3 i — 3/3^ = 0 , 
where from we see that a 3 = 3/3i/( 2 i?o, 2 ) 7 ^ 0 . In this case analogously to the previous case the 
matrix J reduces to the following matrix. 

The columns of the first block are: 

1- st column: ((jo, 20,3 - 6/?i),..., (jo, n -ia n ~ 6/3 m ), 0,..., 0) 

2- nd column: (0, (ji, 3 a 3 - 6/3i),..., - 6/3 m ), 0,..., 0) 


m-th column: (0,..., 0, - 6/3i),..., (j m -i, 2m a: n - 6/3 m )); 

the columns of the second block are: 

1-St Column: (0, . . . , 0, j'o,m+l/^l> jl,m+2^2i • • • 1 jm— 1 ,2 mftrn) 

2- nd column: (0, . . . , 0, jo,mPl, Jl,m+1^2, • • • ,jm-l,2m-lPm, 0) 

3- d column: (0, . . . , 0, jo,m-lPl, h,mh, ■ ■ ■ , jm-l,2m-2Pm, 0, 0) 

(n-2)-th column: (j 0j2 /3 1 , • • • ,jm-i,m+iPm,0, • ■ •, 0), 

where the rn columns of the first block denote derivations with respect to Pi, ■■ ■ ,P m , an d the 
n — 1 columns of the second block denote derivations with respect to a n ,..., ce 3 . 

Now the situation differs from the first main case. If we apply the Gauss elimination algorithm 
to kill all terms of the right part of the matrix lying over terms containing (5 m , we can destroy 
the top non-zero terms. So, we must control the changes of these terms modulo e 2 . Fortunately, 
it is not difficult: the term j m _i_fc in _i_fcO: 3 — 6/?i , where 0 < k < m — 1, will be changed to 
the term 

— 6/3i- ^°'" t+1 k {jm-l-k,2m-kOt n ~ $)Pl = 

Jm— 1,2m —k 

(—1 + m) 2 m (2 + 5m + 2m 2 ) + 2 k 2 (—1 + m 3 ) + k (4 + m — m 3 — 4m 4 ) 

m 3 
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As it can be easily checked, the numerator can be equal to zero only for k > m — 1 . Thus, the 
rank of J is equal to 2m — 1 = dim Y , a contradiction. 

Now we can use the induction: suppose we have proved that the whole component belongs 
to the intersection Y = {cto = ... = cq_i = 0 = (3q = ... = P 1 - 2 } ■ Then the 2(1 — 1) — 1-th 
equation in (1161) implies cx-ifii -1 = 0 . If there is a smooth point with ai 7 ^ 0 , then we can apply 
the arguments from the first main case to show that the matrix J has the maximal rank equal 
to the dimension of Y , a contradiction. If there is a smooth point with /3;_ 1 7 ^ 0, then from 
2(1 — 1 ) -th equation we get 

3 

ai +1 = - Pi -ij 

31-2,2(1-1) 

and, analogously to the case a .2 = 0 , j3\ / 0 , we can control the changes of the top non-zero 
terms (jm-i-k^-i-k+l^+i — 6/3/_i ) , l < m . They will be changed to the terms 


jm—l—k ,m— 1— k+i a i +1 — ®Pi-i 


Jl—2,m—l—k-\-l / . p 

■ \Jm—l—k,2m—k^n ^)Pl —1 

3m— 1,2m —k 


- --r —t (—1 +1 — m) (—4k — 2k 2 +12 kl + 4k 2 1 — I2kl 2 — 2k 2 1 2 + 4kl 3 — 2m + 5km + 2k 2 m + 6lm 

(-1 + 0 "i 3 

— 10 klm — 2 k 2 lm — 6l 2 m + 5 kl 2 m + 2/ 3 m+ 

3 m 2 — 5km 2 — 2 k 2 m 2 — 6lm 2 + 5 klm 2 + 3 l 2 m 2 + 3m 3 + 4km 3 — 3 lm 3 — 2m 4 ). 

The last expression is equal to zero only for k = — 2 + 21 + m > m or 

— 1 + 21 — l 2 + m — lm + 2 m 2 

2(1 — 2 1 + l 2 — m + lm + m 2 ) 

But the last expression can not be integer. Indeed, the great common divisor of m and 
(l — 21 + l 2 — m + lm + m 2 ) must divide also the numerator, i.e. the doubled fraction must 
be integer. On the other hand, it is clear that the fraction is positive and less than one. It also 
easy to check that it can not be equal to 1/2 . 

At the end we obtain that the whole component belongs to the intersection Y = {«o = 

... = a m = 0 = Pq = ... = p m - 1 } with dim Y = 2. Then from (2m — 1)-th equation we 
obtain a m +i = 0, i.e. the component lies in Y 0 {a m+ \ = 0} , whose dimension is one, a 
contradiction. □ 

Let us prove Theorem 1.1 The intersection X' (in P 2m+ ) of X from proposition 13.11 and 
the cubic defined by (1141) is again non-empty, and each its irreducible component has dimension 
greater or equal to 2; the intersection X' D Z with Z is non-empty and each its irreducible 
component has dimension greater or equal to 1. The homogeneous part of (|14D not depending 
on u is 

a 0 p 2 + Pi/ 2. (18) 

It also has a solution of the form P from proposition 13.11 

To prove Theorem 1.1 it is sufficient to prove that for any fixed m > 0 any irreducible 
component of X' (1 Z containing P has dimension 1. 

Note that if ao ^ 0, then either Pi or Pq is not equal to 0. Indeed, if Po = pi = 0, then 
from 0-th equation it follows that p 2 = 0, from 1-st equation it follows that P% = 0 and, by 
iteration, p m = 0 , a contradiction. 

Let Q be a smooth point on some irreducible component of X' f] Z as in the proof of 
proposition ^. II Consider the new Jacobi matrix with the first row consisting of partial derivatives 
of the equation (fl8]) : 

(3$, 20001, 0, ft 2 ). 
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If «o 7^ 0 , then it’s easy to see that this row and all other rows of the old matrix J are linearly 
independent, i.e. the dimension of the component is one. 

If ao = 0, we can literally repeat the arguments from the proof of proposition 13.11 Indeed, 
as we have already seen, in this case even an irreducible component of X n Z would be of 
dimension less or equal to 1. Theorem 1.1 is proved. 


3.2 Proof of Theorems 1.2, 1.3 

According to Theorem 1.1 an arbitrary equation 

T 2 = X 3 + c 2 X 2 + c\X + c 0 , X, Y e Ai 

has infinitely many solutions of the form L^ m = (d 2 + V m (x)) 2 +W m (x), and the equation 

Y 2 = x 2a+l + c\ y X 2g + ■ ■ • + 4x + 4 , X, Y G A 1 

also has infinitely many solutions of the form ip\]Jjf), ^(.Z^ +2 ), where 

L(r) = p\l\) = ((1 - y 2 )dl - 3 yd y + aT r (y ) + bf - ar 2 g(g + 1 )T r (y), 

T r (y) is the Chebyshev polynomial of degree |r| (see section 2). To prove Theorem 1.2 and 
Theorem 1.3 it is enough to prove that at r > 10 and r ^ r\ 


for arbitrary p G Aut(A\) . This facts follow from the following lemma. 

Lemma 3.1. Consider a family of operators of order four with polynomial coefficients 
L(r) = ( a{x)dl + b(x)d x + c r (x)) 2 + d r (x), r G N, 
where a(x),b(x) are polynomials of fixed degree such that 

dega(x) > deg6(x), degc r (x) = r, r > deg d r {x). 

If r > dega(x) + 8, then 


p(L(r)) ± L(n) 


at r / n for arbitrary p G Aut(Ai) . 

Here we assume that deg6(x) = —oo if b(x) = 0. 

Proof. Let us assume that there is p G Aut(Ai) such that at r > dega(x) + 8 we have 
p(L(r)) = L{r\) for some r ^ r\ . Let 


<p(x) = qn(x)d x H-h qo(x), p(d x ) = p m {x)d™ H- \~Po(x), 


where qj,p s are some polynomials. First consider the case n = 0. If n = 0, then m = 1 
otherwise the operator p(L(r)) has order greater than four. Further, 

p(a(x)dl + b{x)d x ) = a(q 0 (x))(pi(x)d x + p 0 (x )) 2 + b(q 0 (x))(p 1 (x)d x + Po{x)) = 

a{qo{x))p\{x)d 2 x + a{qQ{x)){jp 1 {x)p' l {x) + p 0 (x) + b(q 0 (x))p 1 (x))d x + 

+a(qo(x))p 0 (x ) + b(q 0 (x)) + b(q 0 (x))p 0 (x). 
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From our assumption it follows that 

a(qo(x))pl(x) = a(x), a(qo(x))(pi(x)p[(x) + po(®)) + b(q 0 (x))pi(x) = b(x). 

Hence from the first identity we get that p\{x) is a constant and qo(x) is a linear function. 
From the second identity we get that po(x) is a constant, otherwise the degree of the left hand 
side is greater than the degree of the right hand side. Thus 

<p(x) = six + s 2 , p(d x ) = s 3 d x + s 4 , Sj £ C. 

From this we obtain ip{L(r)) ^ L{r\). 

Let us consider the general case n ^ 0. We have the following identities for orders of 
differential operators 

ord <p(a(x)d x ) = ndega(a:) + 2m, ordip(b(x)d x ) = ndegb(x) + rn. ord<^(c r (x)) = rn. 

Let us note that 

ord ip(a(x)d x ) = ord <p(c r (x)), 

for otherwise, since ord <p(a(x)d x ) > ord ip{b{x)d x ) we have 

ord<^(a(x)(9^ + b(x)d x + c r (x)) = ordcp(a(x)d x + c r (x)) = max{rn, n deg a(x) + 2m} > r, 
and therefore ord cp(L(r)) > 2r > 4, a contradiction. Thus, 

ndega(x) + 2m = rn. (19) 

By direct calculations one can check that 

ad(— x) 3 (L(r)) = [[[L{r),x\,x\,x\ = 24 a 2 (x)d x + 12 a(x)b(x) + 12 a{x)a'{x), 


hence 

On the other hand, 


ord<^(ad(—x) 3 (L(r))) = 2ndega(x) + m. 
ip(ad(-x) 3 (L(r))) = ad (-ip(x)f(p(L(r))). 


We have 

ord [ip(L(r)),ip(x)\ < n + 3, ord[[<^(L(r)), p(x)\, <p(x)\ <2n + 2, 

ord [[[tp(L(r)),ip(x)\,ip(x)\,ip(x)\ < 3n + 1. 

Thus, using (fT9ll and our assumption r > dega(x) + 8 , we get 

o Ti 

3n + l > ord[ad(— <p(x)) (<^(L(r)))] = 2ndega(x) + m = n(r + 3dego(x))/2 > — (8 + 4dega(x)) = 

4 n + 2ndega(x). 

We get a contradiction. □ 


Hence Theorems 1.2 and 1.3 are proved. 
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